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Abstract 



We present results of numerical simulations of the kinetics of exciton- 
exciton annihilation of weakly localized one-dimensional Frenkel excitons at 
low temperatures. We find that the kinetics is represented by two well- 
distinguished components: a fast short-time decay and a very slow long-time 
tail. The former arises from excitons that initially reside in states belonging 
to the same localization segment of the chain, while the slow component is 
caused by excitons created on different localization segments. We show that 
the usual bi-molecular theory fails in the description of the behavior found. 
We also present a qualitative analytical explanation of the non-exponential 
behavior observed in both the short- and the long-time decay components. 
Finally, it is shown that our theoretical estimate for the annihilation time of 
the fast component is in good agreement with data obtained from transient 
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absorption experiments on J-aggregates of pseudo-isocyanine. 
PACS number (s): 42.65.Pc, 36.40.Vz 
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I. INTRODUCTION 



Exciton-exciton annihilation is an important process that strongly influences the opti- 
cal and opto-electronic properties of materials at high excitation densities. In particular, 
exciton-exciton annihilation affects the nonlinear and lasing properties of organic systems, 
such as J-aggregates and polymer films.0 Though the importance of this process is well- 
recognized, the microscopic understanding of the annihilation kinetics is still rather poor. 
This holds especially under the conditions of strong exciton derealization and (or) low tem- 
perature, where the usual bi-molecular theory of exciton-exciton annihilation is expected to 
break down. The aim of this paper is to study the annihilation kinetics in weakly disordered 
one-dimensional Frenkel exciton systems, where the exciton coherence size can be consider- 
able (tens of lattice units). This study is of relevance to the optical properties and exciton 
dynamics in J-aggregates and molecular antenna systems. 

The standard approach to describe the kinetics of exciton-exciton annihilation relies on 
the bi-molecular rate equation, in which it is assumed that the effective annihilation rate is 
proportional to the exciton density. This equation reads^i 

h = —•yn — an 2 , (1) 

where n represents the average exciton density, understood here as the number of excitations 
per molecule, 7 is the single-excitation (radiative and non-radiative) relaxation rate, and a 
is the co-called annihilation constant having here the dimension of 1/time. The effective 
rate of exciton-exciton annihilation in the system is then indeed proportional to the average 
exciton density n and is given by an. The solution to Eq. ([3]) reads 

U 7eT* + an (e^ - 1) ' { ' 

where no is the initial population of excitations. If the annihilation dominates the single- 
exciton relaxation (ano ^> 7) the excitation population decreases according to a hyperbolic 
(non-exponential) law 
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The typical picture that one commonly has in mind when modeling the annihilation 
process as is done in Eq. (|1|), is as follows. First, it is usually understood that, as a result of 
strong disorder and (or) high temperature, the Frenkel excitons represent, in fact, molecular 
excitations.lhi Next, it is assumed that the excitations (diffusively) move over the system. If 
the diffusion rate is large compared to the rate of nearest-neighbor annihilation, denoted by 
Wq, two excitations annihilate each other (by fusing into one high- lying molecular excitation 
that quickly loses its energy by vibrational relaxation) when they have reached neighboring 
molecules. In this case a = wq. On the other hand, if wo dominates the diffusion rate, 
the annihilation event may occur at a distance large compared to the the nearest-neighbor 
separation. In this case, the annihilation constant is determined by a convolution of the 
annihilation rate with the pair correlation function of two excitationsjll and may in principle 
depend on time. 

The bi-molecular rate description has limitations, which become rather important at low 
temperatures. Obviously, the bi-molecular approach does not account for the fact that the 
excitons in J-aggregates are generally quite strongly spread, with typical coherence lengths 
of several tens of moleculeslH3 Even at room temperature, this length is of the order of 
10 moleculesjil~0 owing to the large inter-molecular excitation transfer interaction in these 
systems. It is an open and interesting question whether this finite length may be accounted 
for by interpreting a in Eq. (H) as an effective annihilation constant. In this paper, we will 
address this question for low-temperature exciton systems. 

We will study the kinetics of exciton-exciton annihilation of one-dimensional excitons 
that are weakly localized by static disorder. We will use the same framework as was done 
in Refs. [18| , |T9| to calculate annihilation rates. The new element of the present paper is 
to use these rates to follow the kinetics of annihilation. An important step in describing 
the annihilation of extended excitons, is to distinguish between inter- and intra-segment 
annihilation .llHl The rationale for this distinction is as follows. As appears from numerical 



simulations of disordered exciton chains, the exciton states residing close to the bottom of 
the exciton band (the region that dominates the optical response) can be classified into 
groups of a few (two or three) states. The states within each separate group are all localized 
on the same segment of the aggregate, with a typical size N* (often referred to as the number 
of coherently bound molecules), while the segments corresponding to different groups do not 
overlap. In fact, it turns out that the two or three exciton states within each such group 
are very similar in structure and energy to the lowest two or three states that exists on an 
ordered chain of length jV*.&~il In particular, the lowest state of such a group has a wave 
function spread over the segment without nodes and can be interpreted as the local ground 
state. The next higher lying state of the group has a well-defined node and looks like a first 
local excited state, etc. The energy difference of the local ground and first excited states 
agrees well with that of a perfect chain of length N*. 

Obviously, to describe exciton-exciton annihilation, one should consider at least the two- 
exciton states. As is well known, one-dimensional Frenkel excitons are weakly interacting 
fermions (see Refs. P^-f2^). Thus, the wave functions of states with two excitons can be 



composed of Slater determinants of two one-exciton wave functions. Under the condition of 
weak localization, two different types of two-exciton states then appear: (i) those with two 
excitons belonging to the same localization segment, and (ii) those with the two excitons 
localized on different segments. This immediately leads to the distinction of intra-segment 
and inter-segment annihilation as fundamentally different annihilation channel-Hll 

This paper is organized as follows. In Sec. [HJ we present our microscopic model of anni- 
hilation, express the annihilation rate in terms of the basic interactions and wave functions, 
and make the formal step towards the annihilation kinetics at low temperatures, where the 
diffusive motion of excitons towards each other may be neglected. In Sec. |T| we use the 
distinction between inter- and intra-segment annihilation to derive qualitatively analytical 
expressions for the low-temperature annihilation kinetics. A more detailed study is presented 
in Sec. [TV], where we basically exactly solve the kinetics, formally defined in Sec. |TJ, through 
numerical simulations. In Sec. [V] we summarize our findings and discuss the relevance to 



experimental low-temperature annihilation data. 

II. TWO-EXCITON ANNIHILATION MODEL 
A. Motivation 

Under usual experimental conditions, only a small part of the localization segments on 
molecular aggregates are excited. For example, the authors of Ref. 0,[lj§ estimated that 
in their experiments, one hundred molecules per aggregate were produced at the highest 
excitation power applied (0.98GW/cm 2 ). As a physical aggregate normally consists of ~ 
10 4 molecules,!'! while the typical localization segment in their particular case counted 20 
molecules, these authors concluded that less than one exciton was created per segment of 
localization (on average, one exciton per five segments). A simple consideration based on 
the Poisson distribution for the probability of finding an integer number of excitons per 
segment, shows that on the physical aggregate about 80 out of 500 segments are expected to 
be singly excited, while only 8 are doubly excited. Triply (and more) excited segments are 
almost absent. Bearing in mind that excitons, created on the same segment of an aggregate 
or on closely spaced separate segments, will annihilate first, we conclude that a two-exciton 
model of annihilation seems to be quite reasonable as a first step. 

B. Rate of two-exciton annihilation 

As a working model, we adopt a linear chain of M three-level molecules as depicted in 
Fie;. |I|J1B@ The two lower molecular states, denoted "0" and "1", are assumed to form 
an exciton band, as a result of a sufficiently strong resonant dipole-dipole inter-molecular 
coupling. The corresponding exciton Hamiltonian, taken in the nearest-neighbor approxi- 
mation, reads 

N-l N 

H ex = ~UJ2 (KnKn+l + ^ n+ lM + E Elnbtnhn , (4) 
n=l n=l 
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where — U < is the nearest-neighbor hopping integral chosen to be negative, as is the case 
for J-aggregates, and &i~ n (&i n ) denotes the Pauli creation (annihilation) operator of the first 
excited state of molecule n (the state with all the molecules in their ground states serves 
as the vacuum state |0) and has zero energy). The second term in Eq. (^) represents the 
Hamiltonian of non-interacting molecules, in which E\ n = E\ + A n is the energy of the first 
excited state of molecule n with E\ and A n being, respectively, the mean value of the energy 
and a static random offset. The latter simulates on-site (diagonal) disorder and results 
in localization of the excitonic states. Fluctuations of the nearest-neighbor coupling are 
neglected. It will be assumed that A n is distributed uniformly within the interval [—A, A], 
so that the typical magnitude of the disorder is given by the standard deviation o = A/ ^/3. 
For a <C U the exciton eigenfunctions are localized within rather large segments of the chain: 
1 <C N* <C AT-iiHiHiil Throughout this paper, we will assume that this condition holds. 

The high-lying electronic- vibrational molecular term, depicted as "2" in Fig. [I], serves 
as the intermediate state through which annihilation occurs.iSEl We consider one of the 
electronic-vibrational levels to be resonant with the two-exciton optical states and to un- 
dergo an efficient phonon-assisted relaxation to the ground vibronic state (see Fig. [I]). The 
annihilation process itself consists of transferring the energy of two excitons to the high- 
lying molecular term. Assuming this step to occur due to the resonant dipole-dipole inter- 
molecular interaction, we may write the corresponding Hamiltonian as follows 

1 M V 

H *=oJ2 i ^b ln bi m (bt n + bt m ) + h.c. , (5) 

2 \n — mr 



m,n=l 



where V is the matrix element of the annihilation operator for nearest neighbors and &2n(&2n) 
denotes the Pauli creation (annihilation) operator of the high-lying state of molecule n. The 
operator (|5]) annihilates the two excitations occupying molecules m and n and excites one 
of these molecules in the high-lying state. The implication of H a for third-order nonlinear 
optics of J-aggregates has been studied in Ref. |3T]. 

In accordance with our arguments in Sec. |11 Aj, we will assume that not more than two 



excitons are created by the pump per linear chain. Moreover, we will assume that \V\ is small 



compared to the rate, T, of phonon-assisted relaxation in the high-lying molecular state. We 
may then use perturbation theory to calculate the rate of annihilation. Moreover, the back 
process (exciton fission) can then be neglected. The resulting expression for the rate of 
exciton-exciton annihilation starting from the two exciton eigenstate \/j,v) (1 < \i < v < A/") 
is simply given by the "golden rule" 

w^ = -p(E f )J2\(2n\H a \^)\ 2 , (6) 

11 n=l 

where p(Ef) is the density of final states (hereafter replaced by 1/T). Because of the 
fermionic nature of one-dimensional Frenkel excitons, one can compose the two-exciton 
eigenfunctions as Slater determinants of the one-exciton eigenfunctions: 

M AT 

\pv) =EE 1pnrn;vn\±m,ln) , (7a) 

m=l n<m 

1pfj,m;vn ^Pfim^Pun ^Pfin^vm, > (^-0 

where {<p U n} are the eigenfunctions of the one-exciton problem 

AT 

£ flSVm = E v <f vn . (8) 
m=l 

Here, H£™ = (ln\H ex \lm) and E„ is the eigenenergy of the one-exciton state v. Substituting 
Eq. ( ffE| ) into Eq. @ one obtainslil 

2nV 2 M 



w 



E 



(m — n) 



(9) 



m=l Ln=l 

where the prime denotes that n 7^ m. In particular, for a dimer (TV = 2) only one two-exciton 
state exists and its annihilation rate is given by 

12 ^y 2 , ^ 

= = (10) 

In order to arrive at Eq. ([|) we used the fact that Frenkel excitons are non-interacting 
fermions whenever the nearest-neighbor approximation is used for the hopping integrals. 
They become interacting quasi-particles when including the coupling to far neighbors. The 
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importance of the latter can be estimated through the changes which the long-range terms 
produce in the density of exciton states. It is known that in one- dimensional aggregates, 
the long-range dipole-dipole interactions shift the exciton band bottom by approximately 
20% compared to the nearest-neighbor modelJ^H The smallness of this shift suggests the 
corrections due to long-range interactions to be of a perturbative nature. Indeed, for few- 
particle states, Frenkel excitons are weakly interacting (well-defined) fermions, despite the 
long-range coupling. 

Furthermore, as follows from the results of both numerical simulationsi! and theoreti- 
cal estimates^ of the linear optical properties of disordered Frenkel chains, the oscillator 
strengths of the optical transitions near the lower band edge grow by approximately a factor 
of 2.5 due to the long-range dipole-dipole interactions. This results from a larger exten- 
sion of the optically active exciton states in the exact dipole-dipole model compared to the 
nearest-neighbor model (at a fixed disorder strength). In principle, this effect is not of a 
perturbative nature. It can, however, be included into the final formulas for the annihilation 
rates (see below) by rescaling the number of coherently bound molecules iV* . The above ar- 
guments justify the nearest-neighbor framework as a reasonable approach to describe Frenkel 
excitons. 

It is worth stressing, though, that the dipole-dipole interaction in the annihilation chan- 
nel (H a ) can generally not be taken in the nearest-neighbor approximation, because the 
annihilation of two excitons localized on separate localization segments is determined by the 
coupling to far neighbors. 

C. Low-temperature annihilation kinetics 

As in the present paper we are mainly interested in the annihilation kinetics itself, we will 
neglect any other possible channel of population relaxation, such as the radiative transitions 
from two-exciton to one-exciton states (also acting towards lowering the exciton population), 
as well as a possible multi-phonon relaxation from the high-lying term to the one-exciton 
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states (acting, on the contrary, towards raising again the exciton population). To calculate 
the kinetics of the exciton-exciton annihilation, we will assume that excitons are created by 
a resonant laser pulse, that is short compared to the inverse of the J-band width. Under 
these conditions, the initial populations of the two-exciton state \[iv) is proportional to the 
corresponding oscillator strength F^ u given by 



where D = J2 n =i(^tn + &in) is the chain's dipole operator (the chain length is assumed to 
be smaller than the emission wavelength). 

After the initial creation process, excitons may in principle move over the chain, At low 
temperature, however, the possibility to move is very restricted and the optically excited 
localized Frenkel excitons are practically immobile.El The reason is that at low temperature 
(T < width of J-band), an exciton created in one of the local ground states may move 
to an other similar state only when the latter has an energy lower than the former. The 
typical energy offset between the local ground states is of the order of the width of their 
energy distribution (i.e., the width of the J-band). Therefore, after one jump the exciton 
typically resides in the tail of this distribution. The number of states with still lower energy 
then drastically reduces, giving rise to a strong increase of the mean distance to such lower 
energy states. In fact, already after one jump the exciton has a strongly suppressed chance to 
jump further, i.e., such a type of the spatio-energetic diffusion (towards lowering the energy) 
is stopped rapidly and does not yield a sufficient possibility for two excitons to approach 
each other and annihilate. It is worth noting that experiments also indicate the absence of 
such a diffusion, which would manifest itself in a red shift of the exciton emission spectrum 
relative to the absorption spectrum. The experimental data show that such a Stokes shift 
is either absent or has a small magnitude .I'iEl A similar situation occurs in glasses doped 
with rare-earth ions.0 

Following the above arguments, we will assume that two excitons annihilate from the 
positions where they have been created. Then, the time dependence of the population of 




(11) 
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the excited two-exciton states is given by 

P a (t) = 2(j2U,exp(-w^t)\ , (12) 

\ fiu I 

where f^ v = F^ v j Yl^w P^v an d angular brackets denote an average over the disorder real- 
izations. Note that we have normalized the population such that it equals 2 at t — 0. The 
formula (|T2|) will be the basis of our further analysis of the low-temperature annihilation 
kinetics. 



III. QUALITATIVE PICTURE 

Before carrying out numerical simulations, we first provide a qualitative analysis of 
Eq. (0). Following the arguments concerning the nature of the low-energy weakly lo- 
calized one-dimensional states, we separate the summation in Eq. (|T2|) into two parts: 
p a (t) = P™ tra (t) + P™ ter (t). The first part, P* ntra (t), includes all those terms {/iu} where 
the one-exciton states [i and v are localized on the same chain segment (doubly excited seg- 
ments). The second part, P^ nter (t), contains those terms where fi and v reside on different 
segments. The fact that this distinction can only be made for low-energy states is no re- 
striction, as anyhow these states are the ones that dominate the ground state to one-exciton 
and the one-to-two-exciton absorption spectrum. Using the picture of exciton states on a 
chain of effective length N*, one arrives at the intra-segment annihilation rateSil 

<*°=m^w w °- (13) 

where Wq is given by Eq. (|10|) and the factor 57r 6 /18 ~ 270. 

The second term, P^ nter (t), governs the annihilation of two excitons created on different 
localization segments and is characterized by the rateEl 

N* + 1 

< nier = -^^o, (14) 

where R is the distance between the two excited segments. Note that the rate of the inter- 
segment annihilation scales linearly with iV* + 1. As the wave functions of both segments 
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enter the expression for w™ ter , one might intuitively expect a quadratic dependence on 
iV* + 1. However, only one of the two excited segments, namely the one that passes to 
the ground state |0) in the annihilation process, coherently contributes to ■w l a nter giving the 
factor iV* + 1 (so-called superradiant transition). The transition within the other excited 
segment occurs to the high-lying molecular state of each molecule. It is important to note 
that the latter events are summed incoherently, as is evident from Eq. @, thus preventing 
the appearance of an extra power of N* + 1. 

Keeping in mind that R > N* as well as that the numerical factor in Eq. ( |l3l) is fairly 
large, one immediately deduces from Eqs. and fll~4] ) that w™ tra 3> w™ teT provided that 



N* is of the order of or larger than several units, which is the condition we will focus on in the 
simulations. From this, one expects that the kinetics of the exciton-exciton annihilation will 
consist of two distinct parts: a very fast short-time decay, described by P^ ntra (t), and a very 
slow long-time tail for which P™ ter (t) is responsible. In order to estimate the relative weight 
of these two components, we will take into account the fact that the oscillator strengths of 
double excitation of a particular localization segment and excitation of two different segments 
are of the same order. Then, statistical arguments based on the Poisson distribution seem 
to be sufficient for making estimates. If the initial density of excitations equals n (in our 
case, no = 2/M), the probability of finding a typical localization segment (of length iV*) to 
be /c-fold excited reads 

p{k) = (jk^ e - noN * (15) 
k\ 

Therefore, the probabilities of double excitation of a typical localization segment and of 
excitation of two different segments are given by 0.5(n N*) 2 e ~ n ° N * and 2n Q N*e~ n ° N \ 
respectively. Thus, the relative contribution of the shorter component in P a {t) is of the order 
of 0.25noiV*, which drops upon increasing the disorder strength a (or upon a corresponding 
decrease of N*) and grows with increasing the density of excitations n$. 

Now, we turn to a discussion of the character of the annihilation kinetics. We recall that 
in the bi-molecular model, the decay is non-exponential, as a result of the nonlinearity of 
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the driving equation ([!]). We will argue that in our case, the kinetics is also non-exponential, 
which, however, does not result from a nonlinearity. Let us consider first the intra-segment 
channel of annihilation. Here, according to Eq. ([12]), the nonexponentiality is expected 
even in the absence of disorder (i.e., for a regular chain), because the annihilation rate 
w£ u obviously depends on which exciton states are involved in the annihilation process 
(see the next Section). In the presence of disorder, there is an additional source for the 
non-exponential behavior of annihilation through the intra-segment channel. It originates 
from the fluctuations in the sizes of the localization segments, N. Equation (|i~3|) gives the 
typical magnitude of the intra-segment annihilation rate. In reality, iV* in Eq. (|13|) should 
be replaced by a fluctuating value N. Consequently, the sum over the overlapping states in 
Eq. (|T2| ) can be approximately substituted by an average over a distribution of N, G(N). 
One then obtains 

P™ tra (t) « 2 J dNG(N) exp [-w™ tra {N)t] . (16) 

As follows from numerical simulations,0IIEl]'il the standard deviation of N, [/ dNG{N){N — 
iV) 2 ] 1//2 , is of the order of the mean, N — J dNG(N)N } i.e., the distribution G(N) is rather 
broad. Due to this fact, the resulting nonexponentiality is expected to be considerable. The 
numerical simulations presented in Sec. |IV| confirm this picture. 

The origin of the non-exponential behavior of the inter-segment annihilation is two-fold. 
First, the corresponding annihilation rate, as in the previous case, depends on the size 
of the localization segment N (see Eq. (|14|) ). Thus, fluctuations of the latter will affect 
the annihilation kinetics even for a fixed distance between excited segments. However, the 
character of the inter-segment annihilation is determined mostly by the strong dependence of 
the annihilation rate ([TJ|) on the distance between two excited segments, R. The annihilation 
kinetics caused by the fluctuations of R is simply given by an average of the pair kinetics 
2 exp[— w™ ter (R)t] over all realizations of R. In order to obtain an analytical estimate, let 
us assume that the density of excited segments is low (as is in our case), so that R can 
be treated as a continuous stochastic variable. Then the annihilation kinetics is given by 
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an integral similar to that in Eq. (|16D with G(N) replaced by a suitable /^-distribution 
function, G(R). We will adopt a uniform distribution for R, G(R) = N*/J\f, assuming that 
the probability of finding a segment to be excited is equal to the inverse of the number of 
segments in the chain, M /N*. In evaluating the integral, we will extend the integration 
over the entire positive axes, neglecting thus the minimal distance between two adjacent 
segments as well as the finiteness of the chain. Both approximations are justified at a low 
density of excitations. One thus arrives at 

'5\ N* 



■)intra 



it) 



N* 



| o °° dR (l - e -<" te W)] « 2 [l - T Q ?L(N + l^iwot) 1 ^ 



(17) 



where T(x) is the Gamma-function. Equation ( |lTD is correct provided that the second 
term on the right-hand site is less than unity. This holds in a very large time interval, 
determined by the inequality T (5/6) N* (N+l) 1 ^(wot) 1 ^ 6 < Af. It follows from the stretched 
exponential behavior of Eq. (|17]) that further averaging of P^ ntra (t) over the iV-distribution 
will not change the character of the kinetics and results, in fact, in replacing N by iV*. 

To conclude this section, we note that the ("artificial") quantity P a (t) can be simply 
rescaled to the measurable magnitude - the density of excitons n{t) = P a (t)/Af. Introducing 
the initial exciton density no = 2/jV, we obtain 



n mtra (t) w n 



1 

2 



(£) n N*(N* + l) 1 / f W)V6 



(18) 



IV. NUMERICAL SIMULATIONS AND DISCUSSION 

To study the annihilation kinetics in more detail, we have carried out numerical sim- 
ulations for a chain of length M = 200. For such a length, the mean initial density of 
excitations is n = 0.01. We calculated the one-exciton eigenfunctions cp vn by diagonal- 
izing numerically the Frenkel Hamiltonian (|J) for a particular realization of disorder, and 
then composed two-exciton eigenfunctions according to Eq. (|7|). Using further Eq. (||) and 
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Eq. flllD , we computed the rate of annihilation, w£", and the oscillator strength, F^, for any 
two-exciton state \nv). Then Eq. (|T2| ) was used to evaluate the annihilation kinetics. The 
resulting kinetics was obtained by averaging over 20 realizations of disorder. An increase 
of this number did not lead to considerable changes in the calculated curves. As a time 
unit, we used Wq . The results of the simulations for different values of A/U are depicted 
in Figs. |||-[5] by thick solid lines. 

Figure |2] represents the annihilation kinetics for a perfect chain (A = 0), as well as 
a least-square fit by means of an exponential 2exp(— wot/a) (thin solid line), achieved at 
a = 1.63 ■ 10 4 . The fit clearly demonstrates that the calculated curve can not be matched 
by a single exponential. This unambiguously means that not only the states with \i = 1 
and v = 2, having the largest oscillator strengths, contribute to the sum in Eq. flT2|), but 
the other states contribute a comparable amount. The time scale of the kinetics depicted in 
Fig. H qualitatively corresponds to that calculated by using Eq. ([TB]) with N* replaced by 
AT: at M = 200, one obtains w™ tra ~3x 10~ 5 w . 

We also tried to fit the numerical curve in Fig. |^ by means of the bi-molecular equation ([3]) 
taken in the form 2/(1 + bwot). The best fit was achieved at b = 1.07 • 10~ 4 and is plotted 
in Fig. |2| by the dashed line. At first glance, it seems that the latter almost matches 
the numerical data except, maybe, at the initial stage. However, the fitting constant b, 
carrying, in fact, the meaning of the density of excitations (see the discussion presented in 
the Introduction), underestimates the real value Uq = 0.01 by two orders of magnitude. 

Figures [3] and ^ show the numerical results obtained for disordered chains with different 
degrees of disorder, A/U. In Fig. [3|, we plotted the annihilation decay curves in a wide 
time interval, while Fig. f| presents the initial stages of the annihilation process. From 
the numerical results presented in Figs. |3] and |] several conclusions can be deduced. First 
of all, it is clearly seen that the entire kinetics indeed consists of two well-distinguished 
components: a fast short-time decay, becoming faster as the disorder is increased, and a 
very slow long-time tail. The weight of the faster component is smaller than that for the 
slower one and drops upon increasing the disorder strength A/U . 
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It is reasonable to relate these two components to the intra- and inter-segment channels 
of exciton-exciton annihilation, respectively, in accordance with the qualitative picture dis- 
cussed in the previous Section. Indeed, let us use for the typical size N* of a localization 
segment the well-known estimate@®@'0 

N* + 1 = (3vr 2 -H . (19) 

Recall that in our case o = A/y/3. Substituting Eq. © into ([ID, one gets 

< ntra = ^-J^) 2 ^^0.l(^) 2 W0 , (20) 



506 \UJ \U 



Equation (|20|) gives us the disorder scaling of the intra-segment annihilation rate. Accord- 
ingly, we arrive at the following estimates: w™ tra ~4x 10~ 3 u>o, 2 x 10~ 2 u>o an d 6 x 10~ 2 u>o 
respectively for A/U = 0.2, 0.4 and 0.8. Indeed, these numbers qualitatively match the time 
scales of the fast components (see Fig. |). 

We also plotted in Figs. ^ and |] the least square fits of the numerical data by means of the 
function 2— c(wQt) 1 ^ (thin solid line). One observes that at higher degree of disorder (A/U — 
0.8), when the weight of the faster component is smallest, the fitting function fairly well fol- 
lows the numerical curve over almost the entire time interval of decay. The value of the fitting 
constant c = 0.14 — 0.15 is of the same order as the one deduced from the theory, Eq. (|i~7|), 
according to which it must be 2r(5/6)(iV*) 7 / 6 /Vv r = 2r(5/6)(3 v / 37r 2 £//A) 7 / 9 /A/' « 0.24 (The 
discrepancy probably stems from neglecting N* as a minimal separation in Eq. (|l~7D.) This 
unambiguously means that the inter-segment channel dominates the long-time part of the 
annihilation kinetics. It should be especially stressed that the bi-molecular fits, shown in 
Figs. ^ and ^ by the dashed lines, fail absolutely in the description of the numerical data. 

In order to show the character of the decay (exponential or non-exponential) in the 
case of disordered chains, we depicted in Fig. |5] the log-plot of the calculated annihilation 
kinetics for A/U — 0.2. As can be seen, neither of the two components shows an exponential 
behavior. 
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V. SUMMARY AND CONCLUDING REMARKS 



In this paper we have studied the low-temperature kinetics of exciton-exciton annihila- 
tion of weakly localized one- dimensional Frenkel excitons using a two-exciton static model 
(immobile quasi-particles) with diagonal disorder. Our analysis leads to three main conclu- 
sions: 

i - The entire kinetics consists of two well-distinguished components: a very slow long- 
time decay and a much faster short-time drop. The latter component becomes faster with 
higher degree of disorder. The weight of the faster component is much smaller than that of 
the slower one, and decreases with increasing disorder strength. 

ii - Neither of these two components shows an exponential behavior. 

iii - The usual bi-molecular theory fails in the description of the behavior found. 
These findings are well- understood from the existence of two competing options for two 
excitons to annihilate. The slower component is driven by the annihilation of exciton states 
localized on different segments of the chain, while the faster one originates from the an- 
nihilation of doubly excited segments. Fluctuations of distances between two excitons and 
sizes of the localization segments explain the non-exponential nature of the slower and faster 
components, respectively. 

It is worthwhile to estimate the typical rates of both annihilation channels for existing 
J-aggregates. In order to do this, we need information concerning the parameters U, V, 
and T. For J-aggregates U ~ 1000 cm -1 is quite typicalB0@il Less information consists 
concerning the annihilation interaction V, but as we assumed it to be of dipolar origin 
it seems not unreasonable to take a value similar to U. This is, in fact, supported by 
semi-empirical calculations of higher molecular singlet states of pseudo-isocyanine (PIC) 
molecules.B These calculations indicate a molecular Si — > S2 transition that is similar in 
energy and oscillator strength as the So — > Si transition responsible for PIC's well-known J- 
band. We thus take V ~ 1000 cm" 1 . Finally, we will take T ~ 3000 cm -1 , corresponding to 
a vibrational relaxation time in the S2 state of about 10 fs. Using these numbers, we arrive 
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at w™ tra rv 3 x 10 16 iV*~ 3 s _1 . The corresponding estimate for the inter-segment annihilation 
rate taken for adjacent segments (R = N*) reads w™ ter ~ 10 14 iV*~ 5 s _1 . 

At low temperatures, the quantity N* is found to be of the order of several tensi'0Hli. 
Letting iV* = 20, as was reported in Refs. [12] and [TJ], we arrive at w™ tra ~4x 10 12 s _1 and 
w inter ^ 3 x \q7 s -i_ Note that the magnitude of w l ™ ter appears to be even smaller than 
the spontaneous emission rate of a single molecule, which typically is of a few times 10 s s _1 . 
Certainly w l ^ ter is much smaller than the spontaneous emission rate for an exciton state. 
It is to be noted furthermore that raising N* by a factor of 2 will reduce w l J lter by almost 
two orders of magnitude. From the above, an important conclusion can be deduced: the 
inter-segment channel of exciton-exciton annihilation is in fact ineffective at low tempera- 
tures, because the radiative relaxation is much faster. On the contrary, the intra-segment 
annihilation rate is fairly high and should be viewed as the unique way for two weakly lo- 
calized excitons to annihilate at low temperature. However, since this process occurs only 
for doubly excited localization segments, it will affect the entire exciton population only if 
the number of doubly excited segments is high, i.e., at sufficiently high laser intensities. 

We note that the separation into inter-segment and intra-segment annihilation chan- 
nels was in fact concluded from transient absorption experiments on PIC J-aggregates at 
20 K.@0 The part of the observed annihilation kinetics with a decay time of 200 fs may 
indeed be related to intra-segment annihilation, as is clear from our estimate for w™ tra . Our 
estimate for w™ ter shows, however, that it is unlikely that the second component of the 



kinetics reported in Refs. O and 14, with a decay time of 1.5 ps, may indeed be ascribed to 



inter-segment annihilation (see also the discussion in Refs. ^ and |19|) . 

Our findings concerning the ineffectiveness of inter-segment annihilation provide us with 
a way to control the exciton-exciton annihilation at low temperature. Indeed, recall that 
the local ground and first excited states belonging to the same localization segment are 
separated by the energy offset E\ — E* = 3n 2 U/(N* + l) 2 , which is of the order of the J- 
band width.0 For typical J-aggregates, the exciton radiative rate 7, representing the unique 



18 



relaxation constant at low temperatures, is much smaller than this energy mismatch. Hence 
one may get a large number of localization segments to be singly excited by applying a field 
with Rabi frequency smaller than E\ — E*, but larger than 7. At the same time, none 
of the localization segments will be doubly excited. Therefore, under such conditions, a 
fairly large exciton population may be created in J-aggregates, without being affected by 
exciton-exciton annihilation. 
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FIGURES 

FIG. 1. Schematic representation of all interactions contributing to the exciton-exciton anni- 
hilation process. The interaction H ex , indicated by the dashed lines, forms excitonic states in 
the subspace of the molecular states "0" and "1". Excitons annihilate through a high- lying elec- 
tronic- vibrational molecular term "2" (wio ~ OJ21) ■ The first step of the annihilation process results 
from the inter- molecular interaction H a , which induces simultaneous transitions of the molecule m 
to the ground state and molecule n to the high-lying term. The second step results from fast vibra- 
tional relaxation within high-lying electronic- vibrational sub-levels towards the ground vibrational 
state characterized by a rate V 3> H a . 

FIG. 2. Plot of the exciton-exciton annihilation kinetics obtained from numerical simulations 
for a regular linear chain of 200 sites (thick solid line). The least-square fits by means of the 
exponential 2 exp(— w^t/a) with a = 1.63 • 10 4 as well as by the bi-molecular model @, taken 
in the form 2/(1 + bw^t) with b = 1.07 • 10~ 4 are presented by the thin solid and dashed lines, 
respectively. The time unit is chosen to be Wq 1 (see Eq.(|lO|)). 

FIG. 3. Plots of the exciton-exciton annihilation kinetics obtained from numerical simulations 
(thick solid lines) for a linear chain of 200 sites at different values of the degree of disorder A/U. 
Thin solid lines give the least-square fits by means of the function 2 — ci 1 / 6 at c = 0.21 (A/U = 0.2), 
c = 0.17 (A/U = 0.4) and c = 0.14 (A/U = 0.8). Dashed lines give the best fit using the 
bi-molecular model @. The time unit is chosen to be Wq 1 . 

FIG. 4. As Fig. |3| but now focused on the initial stage of the annihilation process. For this time 
interval, the best fits by means of the function 2 — ci 1 ' 6 were achieved at c = 0.21 (A/U = 0.2), 
c = 0.18 (A/U = 0.4) and c = 0.15 (A/U = 0.8). 

FIG. 5. Log-plot of the exciton-exciton annihilation kinetics for a linear chain of 200 sites at 
A/U = 0.2, demonstrating the non-exponential character of the decay. The time unit is chosen to 
be Wq 1 . 
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